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A cational method 1S presented to preq tebe p62 ee 
meerree—-falling objecteminto deep-sea Sediments ye) seems 11 ao 
Peeeven empirical theories. from the field of SoiPaneenani cs. ie, 
Known hydrodynamic phenomena. The impact velocity of the ob- 
ject and the shear strength profile and density of the sedi- 
wear are assumed to be known. The penetration problem was 
fereeved through the use of a computer by equating the work 
feme during penetration to the energy of the object falling 
meroughn air and impacting onto a modeled deep-sea sediment. 
Wemmoo}ects were simple geometric shapes ranging in weicht 
From 500 to more than 1,000 pounds. The impact velocities 
maid £rom zero to twenty feet-per-second. The results are 
Gomoared with full scale tests and recommendations are made 
memercend the method to a water-sediment interface. The 
meee sSuccessfuliy predicts the penetration of objects into 
weak, Saturated, sediments within the accuracy of the state- 
Or—-the-art techniques for measuring the sediment mechanical 
properties. The impact duration time was observed to be 
Meier vely constant and indepencent of object velocity, Shape, 
Peemvecight implying that it may be a unique property of the 


Seeeamic behavior of a sediment type. 
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fie tN RODUCTICH 


A. PURPOSE "Of -LNV Si GALLON 

Various types of objects will be loSt Whwenie saya wena 
the water column and impact on the ocean bottom as man con- 
@unues to move more freely on the surface and through the 
oceans. Such objects can range from small and very expen- 
Peeernstrument packages to fleet submarines.” It is anipeetoane 
Memiaiow GCxactly how far they have penetrated into the bottom 
@eemng Search operations since complete penetration below the 
feieeece Of the sediments may require a search by other than 
fame! OL acoustic methods. 

mae velocity and free fall attitude of an object that has 
been located and photographed on the bottom may be estimated 
Mmeemeecaie penetration distance if the physical and mechanical 
Pe@perties of the sediment are known. The impact velocity 
Meme tude could be potentially useful in reconstrucicing 
the event and also for determining the amount of energy trans~ 
ferred to the sediment and radiated from the point of impact. 

ime Lecovery method and force necessary to break a lost 
object free of the bottom during a salvage operation will 
depend to some degree upon how far the object has penetrated 
into the bottom sediments. 

fie fields of soil mechanics and foundation engineering 
have been concerned mainly with static situations including 
meen prediction of the very slow settlement ehae a StUrucruse 


undergoes as it is erected. The weight is thus applied to 








the soil in progresSively greater amounts. Deep ocean 
bottom installations may be preassembled, floated to their 
Site, and allowed to free-fall to a position on the bottom: 
The free-fall method of implantment is potentially useful 
meremocean Structures since it decouples the Structure ire, 
Eme influence of Surface moticn and tethering lines. Pre= 
@mweeti1On of an accurate initial penetration depth upon con- 
Meme cnererore, becomes a critical design factor. 

An indication of the sediment's mechanical and physical 
meageertlies can be obtained by observing probes of known 
Meeeical configuration as they penetrate the bottom Since 


Mem@etration depth depends on sediment strength. 


B. ECOPE OF ETNVESTICATION 

mace purpose of this investigation was to develop a meenod 
memesedict the penetration depth of objects dropped into @ 
Simulated deep-sea sediment. Experimental data was available 
to test and verify the analytical method that was developed. 
The full scale experiment employed large, heavy objects and 
it was necessary to let the objects fall through air to 
eibow them to reach a velocity at impact equal to the highest 
expected terminal velocity of free-fall through sea water for 
very dense objects. In addition, dropping the objects in 
Geeeecliminated free fall stability problems as the impact 
Peeereude was effectively controlled and the effects of 
angular momentum on the impact were minimized. The sediment 
used in the experiment was of re-worked continental origin 
characterized by low shear strength and high water content 


typical of deep-sea and pelagic Gepesit.- 
9 





The initial boundaries and limitations that were imposed 
on the study were designed to simplify the complexities of 
mae problem to permit the mechanics to be wumger- toed. 
memrectain the ultimate objective of practical appiaeatienee. 
the method at sea. The conditions are summarized below. 

ime Fhe impact velocity at the top of Ene Ssedimenemeou 
Weasein the range from zero to 20 feet-per-second. Zero 
Pemecity at the top of the sediment column involved releasing 
the object at the sediment Surface and its acceleration 
Mmememen the Soft upper layers before coming to fest. The 
Memex limit of 20 feet per second was believed to be a reéal-= 
2Stic terminal velocity in sea water for most objects. 

Z. Free fall was assumed, or gravity was the only down- 
Peeemcorce acting on the object. 

mee Ob}ECE Shapes were limited to Simple geometric. 
meaeepeesctablished hydrodynamic properties from theoretical 
Pee w lations or empirical tests. Configurations designed to 
Memetrace, Such aS a coring tool with a long barrel and large 
maSS on one end were not considered. 

mee Fhe Sediments were taken to be cohesive, with no 
feeter or internal friction. They were also assumed to be 
low strength, less than four pounds per Square inch, and to 
fewer 1igh water contents typical of marine sediments. 

rae scteugd, Gidonot tnelude the GErects sob any tSserele. 
Pera rter initial penetration. 

oe The manimom factors necesSary to predict een erro, 
@epth were considered to be: the impact velocity; geometry 
and mass of the object; and the shear strength profile and 


femserey Oot the sediment. 10 








GS HISTORICAL, BACKGROUND 

The loss of the U.S.S. THRESHER and the  Oeebten sai en 7 « 
fmaechie location of the hull which had penectratedidecee lr siito 
meny Soft sediments served to emphasize the face Ehat Tatele 
was Known about physical phenomena associated with deep-sea 
sediments and the mechanism of penetration. This disaster 
Mmeempoted the original full scale experiments Carried out by 
em@emeaval Civil Engineering Laboratory and used to verify 
mimemmethod developed in this study. The THRESHER disaster 
Peeve pointed out the difficulties involved in Sampling seqi— 
ments at deep ocean water depths and usSing only a few samples 
Maeeweoo in a laboratory to predict as precisely as possible 
Poemsediment's dynamic behavior in-situ. This study is only 
a step toward understanding all the variables and their 
Pee rouship to each other. 

mee following iS a review of the state-of-knowledgemseiece 
the investigation and the technical background. The variables 
magervyed and the development of other penetration equations 
cee SCcussed. 

ie State-of—Knowl]edqge 

The problem of earth penetration by projectiles is a 

ellassical one of terminal ballistics and an immense amount of 
Meee cheoretical and experimental literature exists on the 
Ewoeyect. Empirical attempts to determine the thickness of 
earth embankments needed to protect soldiers were made as 
mey as 1742 [Young, 1969]. In recent years the treo tall 


penetration problem has been studied to provide protection 


abel 








against artillery, to estimate forces acting. Oneane eae ema 
aeing earth impact such aS would be reqmairca >) weeo1 ean] 
meclear power plants that would have to remain amcacteeyring 
ieeentry and landing on earth, and to determine pnysiceal com 
Meepertiecs remotely by an acceleration-time record transmiceed 
from a projectile penetrometer. The stucies have involved 
moemeerelds of engineering, mechanics, fluid mechanics, geology, 
soil mechanics, aerodynamics, wave propagation, aS well as 
Semer disciplines. They were concerned with the problem of 
pmeeoagy falling with high velocity through air and impacting 

on hard continental soils. Bibliographies of prior investi- 
Soeemons are available [Schmid, 1969 and Young, 1969]. Dynamic 
Pameres in the field of soil mechanics on land are still in 
Beemer inkancy. Whether they can be extended to structures 
mmemene penetration problems in the deep oceans probably will 
meme a question for some time into the future for four basic 
reasons. 

First, the velocity of an object free-falling through 
water will be much less than the velocity of a similar object 
mm At very low velocities, from zero to 10 feet per 
meeeme, the inertia of the object may not be important com- 
Bered to other factors. In water the added mass of entrained 
mimeod must be accounted for, but probably would be ignored in 
Skee 

Secorme continental SOLls are quilvesdi1 fe vene Zeew 
deep-sea sediments. Sediments in the ocean basins typically 


mare very Small grain sizes, often clay Sizes of .O003 mm ana 








Smaller. The rate of sedimentation iS measured fractions of 
millimeters per year with very old layered mate urate 
imelatively close to the top of the sediment Column eee 
Sea sediments are completely saturated with interstititial 
pore water often accounting for half their bulk volume. The 
Meme water is not the pure water commonly found on land, but 
the chemically complicated sea water which causes flocculation 
of some sediment particles and the presence of authigenic 
eameaals near the sea-sediment interface. Lower im the sedi-— 
PemeercOlumn, diagenSis and cementation between particles 
wee LOrming a loose fabric often with the porosity being 
mete vely conStant with depth. The effects of the low 2 
fem~mee Centigrade temperature and high 6,000 pounds-—per-— 
See inch pressure encountered at the median ocean depth 
@meeee an entirely different environment than found on the 
@emerments. The degree of reworking by benthic organisms 
mamas be inportant in determining the mechanical character=— 
Memneee of tie upper jayers of the sediments. 

Myra, when measuring = tne mechanical prlopcreieene. 
the sediments, the amount of disturbance tite thie samp ie 
undergoes aS a result of the Sampling process, handling and 
Preeearatlon for testing, temperature and pressure changes, 
Memorgahic growth or decomposition is difficult to precisely 
estimate and may vary from one sampling and testing technique 
to another [Richards, 1961]. 

Fourth, if the sediments in the upper layers are 


mreated aS a suspension of particles in a fluid, which may 


gS 








be a valid approach to sediment dynamics if the water con- 
tent is over 100 percent, then another problem arises. The 
fluid dynamics of non-newtonian materials have been researched 
mote Conduit flows such aS Slurries in pipe@lines? sute jem a. 
marvestigations of external flow about a body |[Pagyashiwanad 


Robertson, 1969], and unsteady flow at a boundary have been 


made. 
fe Prxor Investigations 
Mjeseware Only a few Lecords Of ODjeCeCes ae. 1 ae 
mea lost and recovered. The instances where things have 


been recovered or photographed to document and determine the 
Penetration depth together with adequate data about the 
feeetencal properties of the sediment are virtually non-= 
feeeeeent. The same lack of data exists in reports written 
peeeeeeciiec Lew cases where penetration was investigated. “The 
fw s CYP Cally set forth equations with empirical consctanes 
Smee Oontain no information about the controlling sediment 
Pama nctcers. 

much literature ana empirical 2nformacion Shaya eo 
Peetemor le adriving. The forcing of a long piling into SOEE 
Peemments 15, however, quite different from the free-fall 


MemMeeration problem. 


Be TECHNICAL BACKGROUND 

All of the variables were studied, then the development 
@emeene previous penetration equations was researched betore 
an attempt was made to solve the penetration problem. A 


review of the known variables for which data exists was made 


14 








to develop a method that would be of use at the present 


mime USING availabe cata ang techniques: 


Le inet Va mie Le Sl involved 


Penetration may be a function cf more™thany2O@Nc12 7 


ems variables and may be described by the following eGCaudue.s 


where, 


M 


M, 


C, 


MES ape Gv S IG w, Ss, = 


thy 
(1) 
ey we O pen 


W 
iMpaet VolLoci iy 
MNaSSeGreehe ODby7ect 

GrOsSsS Sectional area of Ehie chiceE 
depth of penetration 

Suctace cCOuchness Of The object 
Obyect Shape 

MiOnG eGLorTO) OL Une Seatment: 

Young's modulus of the sediment 

BO GSCOMeselabtOlwor- Ele Scaminaime 
Cemeicy fOte tile poet ment 

ities OL Dene tiga: Wem CO me eCcle 


Hamer alooic and chem Gal Ceomsesi ti camer 
the sediment 


Mesecostey Ol Ele Sedimenc 
Shear Sereno, Ol eae Secimene 
Peaieaot PLiveOrr ene Secimene 
Samet ctivity Of “Ene Seaimewe 
water content of the sediment 
Omer rOn Of wie SeCdimene 


Sain Size of the Sedrmene 


i 








SITY 


- poSt depositional changes of the vsedqimen. 


p. ] Ceieiiay ore com WwotLer 

i = viscosity of the water 

h = water depth 
a. - adhesion of the sediment to the object 


sucha livst 1s self explanatory and is mee le<— 214. 
Some of the above properties are related to others, and the 
meeameters may thus be reduced in number by a dimenSional 
analysis [Liu, 1969]. Some of the properties, however, are 
meemamenable to a mathematical analysis, for example the 
Meer SiZe distribution and the gqrain's morphology and 
@eeeemcation in the sediment fabric. A phenomological 
meeeach CLO the problem may, therefore, be required [Muacga, 
mae TO put all of the variables into a system of 
SieerOns may be an impossible task, and for ths reason 
meemmerecOl1CCiLON Of penetration into deep-sea Sediments may 
meee SOmewnat Of an art and not susceptible to an Gxact 
peewee on for Some years. 

The adeep-sea sediment particles accumulated over 
mmeerons Of years occur in many types which reflect unique 
feet cional conditions of the earth's history. The 
mechanical properties vary widely, as subsequent to 
deposition the sediments may have undergone secondary 
cementation, diagenesis, recrystallization, and a degree of 
reworking by animals in the top important layers [Smith, 
1969]. Some problems involved in learning more about deep- 


fm=-)sediments are set forth in the following paregrapns- 


16 








The details of the sea-sediment inteémiacewa nome oma 
Sediment layers in the depths of the ocean ereauonl 7.4... 
Ienown the world over. in the most nada OCGan, Che Neen 
Atlantic, Keller [1967], estimates one core has been analyzed 
Moemoecdiment mechanical properties for every 307000 ssquare 
miles of ocean floor. Less data is available in the Pacific. 
Mm@emvariation of the ocean bottom's nature and the high cost 
mmcroout 1,000 dollars per foot of core recovery and analysis 
Wwe limit any great expansion of knowledge. 

Mucheinput tO the design of foundatvons ten land ee 
Pepe On Dpast local experience incorporated into local 
aemreing Codes. On land, tests can be conducted and equa ai Ome 
oewooed that fit special circumstances in C@rtain areas ae 
Pemerocrably less expense than in the oceans. Al1SoO a past 
[ieeweeay CxXists for land areas of interest, whereas Compara— 
tively little has been known of deep-sea sediments until very 
fee times. Although the horizontal area variability of 
deep-sea sediments may be less than that of continental 
Eeemments, to empirically categorize the seventy percent of 
the earth's surface covered by water and anticipate a return 
in less than the number of years that have elapsed since 
Structures have been built on land would be folly. 

fae ne Develooment of Previous Penetration pquamrone 

ph Meyiewmor the development Of Ppenelratlonmeecquaren. 
emphasizes the importance of the various constants that act 
Memmeontrol the penetration depth. The discussion that Tol- 


lows combines the outlines given by Schmid [1969], and 


ey 








Young [1969]. First and most basic from Newton 9S second 
law, 


Bacio Meco 2) 


imemean be shown that penetration depth is determined sr: 
PeEocity iS Known as a function of time. In@@eneral thices. 
Mermecnec Casc for unsteady motion at a boundary. 

Other classical equations are also based on Newton's 
second law where resistance to penetration is a function of 
mi@meastantaneous velocity, V; 


hte Be | De n 
i. a ae + 4- t > 
MX Cy fy (V) ft. (Vv) ae fh Cyan) (Ge) 


weeeely pDOStulates that during the first phase of penetration, 
operon velocity, the motion is governed by resistive forces 
PieePportionali to velocity to some power, and at the final 

Memes Of penetration the resistance is independent of velo- 
eee Li the functions, fo. f., 


MeemetaitsS and the highest power of n assumed to be 2, then 


(tL. are taken to be 


equation (3) simplifies to; 


tei ee ac BV ot yy (4) 


feeteemat high velocities the forces controlling penerration 
ees eOpOrtional to Vv" OF aWalegous to dftagq force 1n7 iia 
awe AL moderate velocities resistance will be proportional 
meme Tirst power of V, analogous to viscous resistance, ane 
meee last stage, proportional to a constant force indepen 
Seme Of velocity. According to how the constants in equation 
(4) are treated, the following classical equations can be 


G@@rived after integrating twice with the initial conditions, 


cS 








Valves =-0, and final conditions, Vv secme ae eee 





mv? 
Robins-Euler z= al (conStant Lesistangmuome ss (5) 
o4 
M V 
Poncelet Za. oy in (l * — =) ‘(resi etiincmeseme (6) 
x a : ae: 
constant and similar 
to ris Gwretrady 
2 
=o Mi ie. 
Resal i pee Cie ote ) (resisting forces (7) 
Cx 


Sime lar seo: 7 ieee 
ViSCOSsaigy sac 
flurzdedrac) 


GepmeerdestiOn relative to the use of the above equations vests 
eeeeiie Constants, a, B, y, Which are determined by the geq— 
feamey sand mass of the object and the type of medium Ilbeing 
penetrated. Schmid [1969] discussed several empirical 
Memeo On equations from the literature which are vari- 
meememc Of the above classical equations. The constants were 
Peeeemrned from high velocity impactS on continental soils. 
Thomason and others [1968] used Pongele 23s) Cais 
Semeecermine the amount of charge required for an eGxplosive 


aeenor system. 


ae ee 
es eg (8) 


II 





where 


p = W = weight of Dio) Cerise 


a: Nemmal L£rontal, area 
a = constant related to the shatter strength 
Of the medium (tor “maid sa =e c 
b = an inertial coefficient for the material 
fore form Laceor, usually ia 








Evaluation of a and b is a problem and unfortunately 
a comparison of predicted and actual penetration was not 
given. However, the penetration depth's dependence on high 
Pemocity was illustrated. 

Young [1969] presented a semi-empirical solution to 
the penetration problem derived from 160 full scale tests on 
@emeinental soils with impact velocities ranging from LOO to 
1,000 feet-per-second. Penetration was treated as three 
phases: (1) impact, (2) moving through the soil, and (3) 
Semeng to rest. The phenomena of entry is not the same as 
femenog through the soil, but the surface effect of entry was 


meme to be negligible after the object had penetrated three 


Peeyediameters. The basic form of Young's equation was: 
fae = 7 ~ ’ 
Z f, (S) f. (A) f., (W) fy (Vo) f. (N) Ga) 
where 
S ~ a Shape parameter of the object 
& ~"“cCross sectional area of the object 


Lee evOLgin, Of ENeROD I ecE 
Vea aL Veloce Ley 


i= Gonstant “dependent Gone soil, peoperume. 
averaged Over penetration distance. 


Two equations were determined for two ranges of velocity, 


@yenwror less than 200 feet per second and one for greater 


than 200 feet per second. The low velocity equation was: 
| -- 2 
2 = 0.53 (N) (8) (yin (1 + 2 x 10°? v4) (10) 
and nomograms were given for casy solution. All Young's data 


fitted within + 20% deviation. Again it should be noted that 
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the soil property constant, N, could vary by a factor from | 
to 50 and "was arbitrarily assigned by those knowledgeable 
mamene field". 
Tf penetration 1s treated as being resisted by aaa 
meeees Only, the following results according to Sehmiawii oo 
MS ay) (ie) 


momen has as its solution 


Ve —Ye 
ee la eC ) (12) 
a 
where 
Ab 
a M (13) 


When t goes to infinity, the maximum penetration becomes: 


V 
= _O 


Za Cras 


a 

but again resolving the viscous parameter of the sediment, u, 
Piemesctablishing a realistic time, t, become difficult. 

Soiis do not behave as isotropic, linearly. elastic 
Meieettals except possibly at very small strains [{D'Appoclonia 
and Lambe, 1970]. Penetration into a material causes the 
feeeeere Limit of the material to be exceeded and theories of 
Peetwercity Cannot be applicd directly to the problem. 
freeones Of elasticity and plasticity applied to problems 
Of deformation assume idealized stress strain laws which are 
Weve! |! suited to soils and more basically ignore the 
important dimension of time, and assume that the state of 
Seer iS instantaneous. 

Various viscoelastic models, Figure 1, with combina- 
tions of Kelvin and Maxwell models, of a sediment can be 
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assembled and solved by Laplace or Fourier transforms, if 


time dependence is linear. 


clastic rebound 


permanent deformation no permonent 
deformation deformation 


aa 


MAAWELL KELVIN 


Figure 1. KELVIN AND MAXWELL VISCOELASTIC MODELS 


Mewevyer the problem of determining the necessary soil param-— 
Peers 1S non-trivial, and to get a unigue solution, Several 
femes Ot Parameters are required [Schmid, 1966]. 

Modern advances in rheology and the continued develop- 
meme, Ot the finite element method with appropriate yielc 
criteria such as developed by Noh and Wilkins in 1964, and 
eeemanea by Schmid [1969], and the analysis of the initial 
Peeerement problem by D'Appolonia and Lambe, [1970], may be 
Useful but again the problem of relating all parameters and 
determining the stress-strain behavior and yield criteria of 
Peeereologically complicated material such as deep-Sea sedi-— 
ments, may remain elusive for the near future. 

im Vepr weak, high water content, decp-ccean sear 
ments, the penetration problem is somewhat Similar to the 
problem of ship slamming. Studies of ship slamming and 
unsteady hydrodynamic impacts at a water surface have mainly 
een concerned with peak pressure, velocity-time relauvone 


ships as a function of the form of a structure as it penetrates 








a fluid boundary and not penetration depth. Extensive 
Peecrature and model test records exist [Ochi ieee. 

A method to @€Stimate the impact veloeime sir enmen- 
eeoenrved penetration depth was developed by Mandi ancmGiveme 
ie2o4)| after an aircraft was lost and crashed into Soft cellays 
Mam@orthern Canada. The method developed and explicity 
solved an equation of motion for a sphere and cylinder 
Me@eerating at high velocities into a Leda clay 


-(m' * M)X = Mg - F - D (159) 


where 


EF esthe Statac strength of the clay aeeuece ren 
CeePCneeracton Gdepuen om, 


D = a drag force -considered avavmctien oe 
PenecEractlon depen, Vereotey sand samme 


Peeeermental tests of high velocity impacts into clay, 
eestor than 800 feet per second) confirmed that the method 
Pmeeeccaicted penetration depth. The over-prediction was 
Sreececoa Since viscous disSipation and elastic wave initi- 
ation were ignored. 

RELer a detailed examination of mathematLicalimetgod- 
Seelabie for analysis of unsteady motions of an object 
Crossing an air-water interface, Moran [1965], concluded that 
Mee@enot them could yield a uniformly valid approximation to 
PeeoOolution of Surface crosSing and that the best hope for 
obtaining a reliable eStimate was with numerical analysis. 

be Known Variables for which Data EX1Ssts 
Since all of the methods reviewed rely totally on the 


ability of the selected variables and constants to correctly 
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parameterize the penetration problem, it Was meéeeeean mae 
review what is known about the variables. The physical 
lmoperties of the object are assumed to be know mt bubreche 
Pemertotion of the sediment strength is subject Women 
mm@eerpretation Of eGmpirical tests or, in Some instances, 
phenomenological observations. Even though deep-sea sedi- 
femmes are very different £rom continental soils, a review oF 
Mmiee nas been learned about the mechanical behavior of fine- 
grained clays on land and what is known about deep-sea sedi- 
wemeao 1s necessary to assemble a workable solution. Also 
eaeworoolems of sampling and testing at sea cannot be 

iene red . 

intormatcron about “che £Ollowing parameters, adenoma. 
Pereete content, void ratio, shear strength, and sensitivity, 
which describe the eh aitie MI be ncuaor of deep-sea sediments, 
fweeerm CX1St in the literature or may be easily determined 
through laboratory tests which have been standarized and 
Peonemay be routinely performed at sea on oceanographic 
vessels. 

The bulk density of deep-sea sediments is determined 
from routine measurement of the weight of a known volume of 
sediment. Typical values range from 74 to 125 pounds per 
Gubic foot with the most common sediments in the range from 
aemmeer O09 pounds per cubic foot [Keller, 1967]. The bulk 
femerties of North Pacific sediments are generally lower 
Meam those of the North Atlantic. 

The VOoud ratio and porosity are a measure Of the 


volume of sediment in a sample compared to the volume of 
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voids, which are filled with sea water in the case of deep- 
Sea sediments. They are also a measure Gheparutele ro saeeee 
point contacts and bonds [Schmid, 1969]. The measurement as 
easily made with a pycnometer that measures that volume of a 
known weight of oven-dried sediment. 

Sensitivity, the ratio of ne Crigimal mower mes. 
paeae Strength to a remolded and disturbed Shear Strength, 
Pmmnoortant to the penetration problem Since the object 
moves through disturbed sediment during penetration. its 
Meeetecance 1S hard to estimate Since most samples are dis= 
[eocd tO some degree. Mitchell and Houston [1969] indicate™ 
miei in the range of Sensitivity from ce S, tyetealmar 
fee>-sea Sediments, the mechanism that causes Sensitivity is 
chemical cementation, which is not well understood by soil 
Meemenists but obviously oe ae 

For a cohesive sediment, the strength 2S a measure ger 
the sediment's ability to resist mechanical deformation. 
Meee masses of cohesive soils generally fail along a Sur- 
face and hence the strength of a coheSive soil is defined as 
shear strength. For purely cohesive sediments undergoing 
bere ederormations, particle to particle interactions whieh 
Pemta cause an “angle of internal friction” and increase in 
Mm@emctrength of the sediment are generally assumed not to 
exist. However, aS more sophisticated testing techniques 
are developed, this may prove to be an erroneous assumption. 

Typical values of shear strengths were given by Keiler 


MmloG7i to range from 925 to 295 pounds-per-square-inch wid 
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the most common range between .5 and 1.5 pounds-per-square- 
inch. The North Atlantic basin sediments haveveliq ei, 
banner strengths than the North Pacific. Une ~reseueeueu 
@ammonates willwsincrease shear strength, bute few carbene 7-4 
emeee round below 12,000 feet. In the Gulf of Mexico) Bryant 
ememothers [1967], found shear strengths of silty clay 
mera tes of the abyssal plains to range from .3 to lil 
Memmas—-per-square-inch with the pelagic sediments of the 
S@oeanental slopes having strengths up to 5.9 pounds-per-— 
S@igewee—inch. 

Rechards [1967] predicted that axreducti ongoee sim—c1eu 
strength would be caused by sampling and laboratory testing 
Peemenat the reduction might be proportional to the grain 
Seer che amount of gases in solution, and temperature in- 
Sfeaces that cause increased waereriad actizon. Hoewevews 
Seimees Lestcd in-situ at 300 meters with a vane device 
memim@etated well with tests from cores taken from the same 
pie a year earlier and tested with a fall cone and laboratory 
mm~emereachards, 1967}. Recent experience indicates this may 
Meeeeec)l) an isolated instance of agreement. 

iaboratory Cests to determine the ericcts of tenpeea. 
Meme and pressure on shear strength and consolidation have 
been carried out with greatly disturbed, remolded samples 
[Vey and Nelson, 1966]. The results were somewhat incon- 
@imsaive due to the small number of samples and difficulties 
iMm~istrumentation. The results indicated that there was a 


decrease in direct shear strength with increased environmental 


26 








pressure for a fine grained sediment With high wo veer 
ioe vane shear testS with increased environmental  ereccun 


pa@wed an “rcrease in shear strength for plastic Ssedimences 


ence a dec; <NCLASSIP TED ppp "cee nlastic sediments. 
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Siexela problem are 


@imiae. . ‘S and must often 


be physically contained during... ting. This pre- 


Sen - 


Swmees tne use Of an unconfined compression test or any test 
that requires trimming the sample before testing. The vane 
@@ear and fall cone tests can be done on a Sample in its 
Original container and the vane shear has the added advan- 
tage of easy adaptation to in-situ testing. Kravitz [1970} 
compared the results of Wykham-Farrance vane, Torvane, and 


fall cone tests on weak sediments and found that the results 








from the fall cone were more repeatable than them aon. 
barrance vane or the Torvane which showed least tepeaca ae. - 
mmese results may have been expected Since an impose as 

mictOr in vane shear testing 1S the rate of Toeadingmencdmen- 
Peearation Of Static and dynamic resistance [Housel o)|- 

The vane Gest s advantages are: ~Paest > te yeamme- 
Gonducted in-situ and a strength-depth profile is easily 
Seearned where all other methods require reentry of a hole, 
mmm may be very difficult in the deep ocean, or discrete 
Pemeling which disturbs the sediments. Second, it if a 
Simple inexpensive technique when compared to sampling and 
Mmiegieatory testing. Third, it does not require a system of 
emotes or free-falling objects to cause the soil to fail 
Piemenerefore has the potential to be used at sea on a 
meme, Ditching ship. . 

The vane test limitations can be Summarized as: 
(ieee it Can only be used in uniform, saturated, cohesive 
Peemment. Second, when done in-Situ it does not obtain a 
emer for evaluation. Third, it imposes a failure surface 
which may not be the weakest surface or relevant to the 
Eeoowemn being studied. 

Eden [1966] found that the vane will yield higher 
strengths with more consistent results in sensitive clay 
with shear strengths less than ten pounds per Square inch 
fien compared to a fall cone. This is in conflict with the 
more recent tests of Kravitz and emphasizes the problems in 


testing. Kravitz made statistically significant tests ona 








laboratory prepared montmorillonite sample and Eden tested 
Matural clays. 

Shear strength values from a vane were found to be 
15% higher than those obtained by an unconfined compression 
Gest on Gulf of Mexico Slope sediments [Morelock, 1969]. 
Mees may be partially balanced by the findings of Ind@ebitzen 
[1970] who found that laboratory vane measurements were 227% 
a ieechnan in-Situ measurements made from DEEP QUEST at a 
@eotny of 3,500 feet. This is in rough agreement with the 
estimation made by Crisp [1968], of a decrease in shear 
strength due to an expansion of pore water caused by tempera- 
MePSemecic) pressure changes. Crisp estimated that an expansion 
peempere water by two percent would occur by raising a sample 
meom 22,000 feet deep and 2 degree centigrade temperature to 
Peace temperature and ee THIS Gxpanslonmeouta 
fmrmimekcantly alter the strength properties of a sediment 
Bema Crisp Showed shear strength Should decrease 10 to 15 
Bercent. 

The in-situ tests conducted by Inderbitzen showed a 
Mig@eonm increase in strength with depth in the sediment, but 
the laboratory tests were erratic, suggesting fase ee di fies. 
ence may have been caused by disturbance during sampling and 
Peemeeration for testing. 

Penetromerers Offer the great advantageyor porentiaus, 
being self-contained instruments. They either record or 
transmit a force versus time curve from which depth of 


eemetration by integration and soil properties Can be 
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inferred. In order to extrapolate to difiteventeenasear 
velocities, however, the curve's shape must be a character- 
ioietc Of the Sediment parameters. There are reqivereer 
With instrumentation and directional stabilization as the 
Pemietcrometer penetrates since centric impact rarely ween 
Peeeneother than symmetrical objects [Schmid, 1966]. Three 
Sener problems also complicate the penetrometer approach to 
Semetmmining Sediment properties in the oceans. First; to 
wempacep information, the penetrometer must have enough 
energy to penetrate, and at low terminal velocities in water 
meee probably require a Smooth, Known, driving force. 
second, viscous forces may predominate on Small penetrometers 
fees extrapolation to larger objects very difficult. Third, 
Pe wereroneters do not retrieve samples for analysis of other 
Sediment parameters. . 

PeEOUSTIC Sounding and seismic methods Go not discern. 
inate enough of the structure of the upper layers of sedi- 
tem@eamimportant to the penetration problem to allow their use 
RePeemocesent time. In the future the use of high £requene, 
Peeeetee sources mounted on deep vehicles, hold promise for 
Semect resolution of the upper layers. 

Une mechods Of G€etTermi ning Sediment, streng en eeliee 
Memwemocen discussed were developed for use in classical 
foundation engineering where a safety factor based on 
engineering judgement would be applied to the final result. 
MmemocnetLration problem is quite different as it requires an 


Poe answer Since an error factor of two could mean the 








difference between an object being completely buried below 
the surface or having one half of the object exposed and 


easily located. 


Sal 


Lie METHOD LOR ANAS Sie 


Prior to a discussion of the method used to solve the 


penetration problem in the air-to-sediment case and the 


modifications required for use at sea, it iS Suitable to 


quote directly from Terzaghi and Peck [1948]. 


~Because ot the unavoidable undertaintres on ole 
in the fundamental assumptions of theories and in the 
Mierecacal Values of the Soil CoOnStantsS, Simplicity ac 
mach Greater importance than accuracy. If a theery, 
momsample, one can readily judge the practical con- 
Peewence Of Various conceivable deviations from the 
PEinoeLOnsS and act accordingly. Li aw tneory 2s 4e0r 
iereatced, it serves no practical use until the results 
meemcondensed into graphs or tables that permit rapid 
meeiiation Of the fanal equations on the basis of 


Several different assumptions.” 


Pee osUMPTIONS THAT PERTAIN TO THE SEDIMENT 
Deep-sea sediments are not well defined in the literature 
Sealing With soil mechanics. Therefore a list of important 


PIwinocions considered in developing the penetration pre- 


Mme ta1on method are stated: 


1. The sediment is completely saturated with sea water. 


2. The water content of the sediment is 75 percent or 


Meme, es 1S typical of marine deposits in the deep ocean 


baSins. 
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3. Organic material that would tend to bing stiee2-717- ae 
moecticles together is either absent Or 2S presemmminm ea 
Smatil quantities. 

4. The sediment will fail along its weakest plane. 

5. The only shear resistance is cohesion. 

6. The initial penetration takes place in such a short 
meatemenat the effect of permeability is nil. In a Study of 
jame—dependent deformation of clays Barden [1969] noted that 
tome effect of adsorbed water on individual particles was that 
meee pressure equals the applied pressure at the first 
meplication of load. 

fmeenny laithification or cementation That may have Mevetoped 
meen naividual particle-to-particle contactS is broken at the 
mien of impact. This implies a portion of the sediment is 
memOolided at impact. if the partic do not react witiewme 
another by other than mechanical means, the mixture will 
behave aS a viscous fluid [Peck and others, 1953, Pazwash and 
mobertson, 1969]. 

8. Shear strength varies in an irregular manner with 
Meo ene variation of velocity gradients is greater in Ene 
deep ocean basins than over shallow areas [Ewing and Nafe, 
1963) which implies that deep ocean sediments have been 
deposited in layers. This is confirmed by visual observation 
of deep-sea core samples which often show horizontal layers 
@eecittferent colors and textures. When studying in-Situ 
Merengths, Inderbitzen [1970] found that horizontal vari- 


merons were much less than vertical variations. He also 
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found that the percentage of variation in? sheave mwe meen 
Gecreased with depth, which may indicate that biological and 
ememical processes are taking place in the upper ayer ene 
feorimal consolidation occurs. 

9. The bulk wet density of the sediment can be considered 
@emstant with depth in the surficial sediments. The shear 
eet Of a Saturated sediment devends to a Targe extent ‘on 
fme Structural arrangement of the particles and Scott [1967] 
Poermmes out that soils can have different Structures at the 
Segemoensity, but may not have the same structure at diftfer— 
ent densities, which implies that shear strength is not a 
mere hunction of density. X-ray diffraction and polarized 
Meee microscope techniques cannot reveal the details ©f 
mame le configuration or orientation in fine-grained marine 
Pea@ments. However, electron microphotographs of Seaimemre 
from the Gulf of Mexico Show them characterized by a loose, 
Been arrangement of randomly oriented particles and that 
MeetrciG rearrangement resulting from gravitational consoli- 
Meleron 1S not apparent [Bowles, 1969]. 

ire The energy transmitted to the sediments ihn the form 
Mewelastic or gravity waves from the initial impact shock is 
negligible in comparison to the total energy dissipated during 
penetration. Impacts on soils and fluid media have been 
meiared for a long time, but there is a very difficult instru- 
mentation vroblem with detecting energy transmitted away from 
Siesimpact site. According to Hunter [in Schmid, 1966] the 


Hertz impact theory, which refers to impacts in the elastic 
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mange, States that the portion of total emercy tia 
Paleieational energy is very small as long aGeyene impact 
Temrocity is less than the velocity of sound in the mecivun. 
Schmid [1966] meaSured the diStant dynamic disturbances 
caused by a small penetrometer impacting into materials whose 
mieological properties were known or easily, tested, such as, 
Meeeetees, Wax, and homogeneous soils. High frequency vibration 
meso ecetected in the hard materials and it was determined by 
MeeeOoL Surface accelerometers that the largest amount of 
energy propagated from the point of impact was 4.3 percent 
@mmmeiie cCOtal energy available. His calculated propagation 
meeea suggested that the motions measured at the surface 

were Stoneley waves. Ochi and Bledsoe [1962] when conducting 
Mmeagel tests of ship hull forms slamming onto a water Surface 
Pounce that forces were very large for .0Ol1 seconds and presemee 
MeewtLO impact phenomena were over within .Q7 seconds, yet the 
model's velocity decrease during this time was less than 4 
Meeeene. Chuang [1970] in a similar investigation, did not 
Seeeet any Significant acoustic pressure variations at a 
distance away from the impact point and concluded that water 
Mould be considered incompressible for impact velocities from 


Metro 50 feet per second. 


B. METHOD OF SOLUTION FOR THE AIR-TO-MUD CASE 

A rational, phenomenological approach was used to solve 
the penetration problem rather than an empirical method 
because of the poor state of knowledge that exists relative 


to actual behavior of deep-sea sediments. With engineering 
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judgment, this method may be extended to different objects 
ema different Cases more easily than an Gmpirical) aopee-. nese 
moe. problem. 
ito Role of Hydrodynamics 

The sediment was treated as a homogeneous, income eese i. ae 
fluid to enable an estimate to be made of (1) the momentum 
MmeamoeLrerred from the body to the sediment, and (2) the drag 
forces exerted on the object at high velocities. If the sedi- 
ment was a true fluid it would have an infinite Strain for an 
Peerred Stress. Therefore, the forces required to initiate 
failure in the sediment and bring the body to rest were taken 
Beecne forces calculated using the ultimate bearing capacity 
Pormula Ceveloped by Prandtl and Terzaghi. 

mrCataumeg Gee@e-Sea sediments as flmids has beenean 
Pepeoach applied to the study bf turbidity fELlows anc sSlumporme 
Poecie COntinental slopes. The term “Spontaneous liguefaction 
Mmenoften used when referring to saturated clays and implies 
mera properties. Briefly, Morgenstern [1967] states that 
Mrevetaction occurs when the sediment is loaded to the collapse 
Stress of the sediment fabric and the load is then transferred 
memene interstitial pore water. The increase in pore pressure 
then produces hydraulic gradients which further disrupt the 
sediment fabric causing the solids and pore water to flow 
Meleroscopically as a viscous fluid. Houston and Mitchell 
[1969] describe the fluid phenomena as occurring in sensitive 
meio which contract Guring shear because of their open 


meric. If the clays are saturated, the contraction tendency 
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meesults ina large transfer of normal Stress from toewele) 
skeleton to the pore water cauSing pore pressure to increase. 
In saturated deep-Sea sediments this transfer of normal stress 
memos almost inStantaneous. Continuing the descr iperoaees 
Houston and Mitchell [1969], as the remoulding caused by con- 
menuing Strain becomes more pronounced, the effective normal 
@me@eoo may fall to a very low value. Very little Strain is 
meemired to break bonds at particle contacts and there is 
m_eeemomi|O-Up of pore pressure even before a significant 
memoer Of bonds are broken because of slight elastic compression 
Seems fLabric and incompressibility of the pore water. The 
merstance of sensitive clays due to dilational energy was 
Mettermined to be small or absent by Houston and Mitchell [1969], 
and Walker [1969]. Because, first, at high void ratios little 
or no tendency toward Pune eee ae observed Quring Shear even 
Memo strains, and second, the effective stress had been 50 
meeneed at high strains the mobilized resistance of the clay 
fabric was small. Mandl and Givens [1964] when conducting 
Meron velocity penetration tests into a clay, judged that the 
Compressibility phase of impact lasted until the surface 
Boundary breaks and that the static yield strength of the 
clay was exceeded by many orders of magnitude under impact 
eng the clays behaved in a dynamic manner Similar to a fluid. 
Schmid [1966] when conducting low velocity free-fall IMpaAc st 
experiments on a clay also observed that liquefaction oceby rear 
Another reason to treat deep-sea sediments as a fluid 


was that their bulk wet density has a narrow range of variation 
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with depth in the upper layers of sediment and the assumption 
©£f a median density over the range of expected penetration is 
a good approximation. 
mee Ultimate Bearing Capacity Equatiog 

The bearing capacity of a Soil 1S eaa@imeasurere ass 
Baeacy CO Support a load relative to its limiting ecuilibeium 
Mmeeording to Harr [1966], the first rational approach to 
Bere Oping a method to predict ultimate load carrying capacity 
See sOll was provided by Terzaghi who recognized that if the 
ees dominant factors contributing to bearing capacity, the 
Gameiete OL the soil, the effects of surcharge, and the Strength 
Parameter of the soil, could be assessed, a good estimate of 
Peeing Capacity would be obtained. Terzaghi chose as a model 
Oe the failure mechanism of the soil the solution to the prob- 
mem Of penetration of Paco: punemes into metal as developed 
PeyePrandtl and Reissner. In the original solution Prandtl 
Was able to neglect the effect of the weight of the metal 
meveaq Co the influence of surface loadings on the surrounding 
material. When the material was assumed to be isotropic and 
mromogeneous with no angle of internal friction, Prandtl’s 
Selutron yielded, 


dg = yD, = 5. Las (16) 


Oa the basis of experiments and Prandtl's solution, Terzaghi 
and Peck [1948] derived the following semi-empirical equation 


Mmerethe bearing capacity of a circular footing. 


Va ae PD: N, ot eye (7) 








where q —~ pressure at Surface to cause failure | oca) 
@Gs> cohesion or shear strength of Siem somi oo 
ee) eee ete he Ot Ve iecumeonel 
D. = depth of footing below the surface of the soil (in) 


factors which depend on the ane@e ror internal 
Perec tloumceuic: CODthem HoOOe iMac 
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The first term represents a contribution due to the 
Sear Strength of the soil, the second the effect of material 
removed by putting the footing below the surface (equivalent 
to buoyancy in water), and the third term, the effect of an 
angle of internal friction (assumed to be zero for deep-sea 
sediments). The equation is modified slightly for Square or 
Oblong footings and Meyerhof [1961] gives constants for wedge 
and cone shaped footings. More complete discussions of the 
Bearing capacity equation and slight variations in the values 
Seeeelye COnStants can be found in texts on Soil mechanics 
feeZaghi and Peck, 1943; Harr, 1966; and Hough, 1969]. 
Terzaghi's equation was the closest to Prandtl's (which 
requires theoretical failure Surfaces) and it has been 
@m@encrcally verified for many years. It was therefore 
Selected as a reliable relationship for estimating the static 
strength of deep-sea sediments. The values of the constants 
and the shapes used in equation (17) are given in Table [I. 

lite bearing Capacity equation, “ehe Shear mest ee 
ance mobilized on the failure surface of any Surcharge 
material is neglected, hence the equation is conservative 


ena relates only to events occurring below the footing or 
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penetrating object. buring penetration, static resistance 
will occur in the undisturbed portion of the sediment below 
mene deepest penetrating part, the other surfaces Of vine 
penetrating object are subjected to remolded sediment or 
meaqiment already set in motion by the penetrating oebijece. 
Bmeretore the bearing capacity equation constants used in 
Pemetration prediction were considered to be for shallow 
meerings for all depths of penetration. 

fen new eharacter of Failure 

Before a final method of solution was selected, 
Moreton pictures of penetrations into soft sediments and still 
Bmovcographs of soil failures were studied. In the motion 
Mectures, the sediment was observed to move laterally outward 
from the object and have somewhat of a separation point from 
Bae Object. While dropping plates vertically i1nveraseine 
lated sediment, Erchul [1968] found that polymer coatings to 
meee Orag on the plates had little effect at high entrance 
melocities and postulated that the sediment may not have been 
memmecneact With the plates. He also observed that at low 
Mmeereics the lubricants did reduce adhesion to the plate 
and considerably greater penetration occurred with lubricated 
ertees COMpared to unlubricated ones. 

Tosti qualitative indication of ue Clbeceuna: 
fluid in interstitial spaces, Rowe and others [1962] dropped 
mreobes into finely ground (less than 40 microns) silica dust. 
At atmospheric pressure the dust displaces laterally and 
Sraters were formed. In a vacuum with no air in the inter- 


stitial spaces, resistance to penetration was 100 percent 
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greater and the dust did not move laterally and no mounds 
Sencraters were formed: 

The motion pictures of penetrating Objects in the 
More pit revealed little splatter and 11ftlegsi «at .cegen 
mupeeance at a distance from the penetration Ppormme Oo low« 
me@meRopertson [1965] observed motion pictures of a buoyant 
Peereecxicing from water and came to the concluSion that 
mregy LOSt due to Surface effects was small. Moran [1965] 
@lse arrived at the same conclusion. 

4. Eguation of Motion 

Since penetration iS a dynamic phenomena involving 
Baemunsteady motion of a body, the first attempt at evalu- 
feeeeomcie penetration depth involved writing an equation of 
femmermeror the object as it penetrated the sediment. The 
mepatecOlution required a ere energy appreach thatvallevwed 
Mememeange in velocity and penetration depth to be determined 
SeemeimenLally with a computerized numerical scheme without 
Mma@emeecuirement to integrate the equation of motion and 
eetadolish ae IMS eateaneOus OC lationship leEwCemm mo mee se llaiae 
and acceleration. 

The forces acting on an object penetrating a sediment 
Biearaccerized as a fluid with a static strength are shown in 
fee 2. The equation of motion can be written as follows; 

Mag = > Applied Forces (18) 
SeG@e tie applied forces Bea weight of the cbjece, Ene Scdre 
Ment Static strength, the buoyancy provided by the sediment, 
Eme force to overcome the inertia of the sediment, and the 


arag forces. 
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Figure 2. FORCES ACTING ON A PENETRATING OBJECT 


mnie change of the body's motion as it penetrertcemen— 
Memirent 1S the Gquivalent of the resistance encountered by 
MmiemOOody, or the energy involved in the change is eguai £0 
meme expended in propelling the body. Momentum is trans- 
/Mered to a Surrounding fluid by; (1) a certain volume of 
fara being accelerated uniformly, (2) the fluid being put 
into turbulent motion, (3) compression waves, (4) gravit 


waves or spray [Hoerner, 1958]. All of the momentum is 
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eventually consumed by viscous friction. The Vast e cosmo. 
of momentum transfer have been assumed to be negligible for 
meme velocity ranges and object sizes under consideratrcniecd. 
were discussed previously. 

The treatment of the drag foree, or the Momentum. ms 
maeO Che fluid by turbulent motion and viscous drag, outlined 
below, follows that given by Mandl and Givens [1964]. 

The drag forces cannot be evaluated easily for unsteady 
feeon through a fluid boundary. The coefficient of Grag can 


ie approximated by; 


B 
Maar 
ne Dy 


Ca 7 ( 


ie (1:99) 
Meret irst term represents the contribution due to viscous 
Grag; the second is that due to form drag; 86 is a constant 
depending on the shape of the body and is usually represented 
by the first term in an infinite series (Oseen'sS approximation). 
64 
T 


mea Circular disk placed broadside to a stream £7 a ele ee 


for a Sphere 6 - 24 (Lamb). The term C, corresponds to the 
Pee cCOLIC case aS Reynolds number approaches infinity in 
Navier-Stokes flows, it depends on geometry and not the fluid 
Gynamic properties. Inspection of experimentally obtained 
Grag coefficients versus Reynolds number relationships for 
Smite bodies show that Ca is a Slowly varying function at 
large Reynolds numbers. For bluff bodies with sharp edges 
Moly all the drag is form drag. The experimentally Geter- 
Paned values of Se the case of bodies with Sharp edges isc 
approximently 1.0. For bluff bodies without sharp edges, 


Pien as spheres, the form drag also depends on the position 
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ef flow separation which in turn depends Gh the era 





mange Of Reynolds numbers. The equation ef iiererenma1@ 
BecOmes : 
— 1, Pe 1 2 
t+ = m w ee aye 
(m'*M)x Mg Po ¥ AX 5 A Dek 3 PAC, x (20) 


Lf the viscosity of clay is taken €o be ometierenean 
ee 3 xX ions lbs - ee a [Whitman in Mandl and Givens, 1964], 
meen Reynolds number for the objects that were under study with 
weeearee Loot maximum dimension and with a velocity of ten 
Meeeeoct-SCCONd, iS On the order of 2 x One Admittedly a 
Bem@etant Value of viscosity for a rheologically complicated 
mmoetance Such as a deep-sea clay 1S a gross ovér Simplifi-— 
Seton, but the order of magnitude of the terms indicates that 
Meme sSCOUS Crag contribution to the equation of motion may 
memeromored {Mandl and Givens, 1964]. Also, towards the 
meena! phase of motion where the effect of viscosity is 
ieaole to increase rapidly with decrease in the rate of 
meen the influence of viscosity will be negligible Since 
faewvelocities are small and the other terms in the equation 
are larger by several orders of magnitude.: Since form drag 
represents vressure forces on the body vanes are Signi fea. 
Meme aigher velocities, it could not be neglected. 

wee added Mass Ettects 

Another reason to treat the sediment as a fluid was 
moraave a method available to estimate the volume and velo- 
Sity of sediment displaced as the object penetrated. Womack 


[1967] verified that the amount of saturated clay moved 


during penetration was equal to the size of hole created. 
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The energy required to move the sediment is conserved since 
ence the sediment iS put in motion by the object it will 
moencinue to move away from the object créatang asta. = 
meteer as the object is Slowed by other forces syle mcm oa 
mmeortant though, because the drag term of the equation GF 
Motion is dependent on velocity Squared and the inertial 
ee-ct Of the sediment motion will alter the velocity-time 
MietOry Of penetration which hasS Significant effects depend— 
memememDOn the shape and impact velocity of the object: 

Added@imeass iSeeaused by hydrodynamieworessure Forces 
Meme are proportional to the instantaneous values of accel= 
Peeeron and according to Wendel [16950] Vase eC aise © Oma Ve seine 
same value in a viscous or inviscid fluid as long as small 
moares are not involved. The added masses were qetermincd 
ei@eeording to equations and methods given by Saunders [{195/7], 
femme are tabulated in Table II. 

The assumptions pertaining to the fluid approxima wen 
Saececp-Sea Sediments can be summarized as follows and arc 
m@emscame aS given by Moran [1965] in a Study of water entry 
mae) @xit’. 

ime trrOorebional t£low — viscous effects On -ehewpre. 
sure felt by the body are generally small and the variations 
Mmaeimviscid forces are more important than viscous drag. 

oe) dncompressible flow - impact velloei tree ware vel 
below the speed of Sound. 


Seo ede@e erteecrs are negligible, 
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The problem with determining Gépth Of penectecaeue: 
from equation (20) is that it cannot be integrated during 
mae unsteady motion of an object transiting a fluid boundary 
mae one fluid to another. This is true for the folle imag 
reasons: (1) m', the added mass and mass of the displaced 
Beeament, iS a function of depth of penetration for an object 
whose submerged portion changes shape as it submerges; (2) 
mremconinant term, that of static soil strength, varies with 
m@eraepth in a highly irregular manner; (3) the location of 
Pie air—-to-sediment boundary condition is unknown a priori; 
and finally, (4) the flow is unsteady. 

oe Transformation of the Equation to Work—Energy 

Me equation of mMotLon may easily be resclved@ine.e 
Maeecuation of energy and work which is time independent and 
Peeemnetion of increments of distance. The development is 


eeevern below. 








ce 
ee ae, (21) 
Cie 
Peeeeeolacation by dx gives 
zm q Gxt . mes ax 
> F. dx = M ax Breas Bierce) (22) 
Mimemtirst integral with respect to displacement is: 
r ax = 1am @%) +e ce 
Se : GG 


The left hand term represents the work done on M Gur diner wre 
interval ax and can be evaluated if BP is a known function 
of x. The right hand term represents the corresponding 


change in kinetic energy of the object. 
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Kinetic energy 1S a Scalar quantity that seepenaam. ye. 
@n the mass and magnitude of the velocity. (‘Sinee tne else, 
memcgquared 1t 1S always positive. The topalemechoantea een... 
Seecne body is equal to; 
E> Kear PE (24) 
Memeo Object penetrating the sediment and Coming to rest, 
meewenergy of the object can be accounted for by: 


+ 4 = 4 _ 
1) a our eae (25) 


mes Simple balance states that the total energy E, of the 


i 
Meyect at the beginning of an interval plus any energy put 
[eerie ObJect minus any energy taken out must equal the 


mal energy, E of the object at the end of the interval. 


2! 
When the object has lost its energy in overcoming buoyancy, 
meter sOil strength, inertia of the sediment, and drag, it 
wl Scop its motion. 

Pespartiecular reference level for "the eCnergrec toa 
Mereeo since energy is relative. The potential cnergy 
Peweeros On an arbitrary selection of datum for zero potential. 
Beoo tne kinetic energy is arbitrarily expressed relative to 
Meme aniition at zero velocity. Hence any physical meéasure-— 
Meme Over an interval describes a change in energy. it 
follows that reference to an arbitrary datum used for the 
Seoression of energy will cancel. 

The work-energy method approach to solving a proplem 
is powerful but dangerous because an answer is always obtained 


mre no indication that the answer is correct. The metiice 


balances work to egual energy available and if a portion of 
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ef: the energy available is not cConsidereqs om amcememeemin, 
Benouted, work iS Still done and an answer results. In the 
Mmsteady case of penetration, uSing an iterative solution. 
Peaost requires a computer because errors in compucae: one et 
not be indicated but will appear as energy gained or lost. 
The answer will not appear wrong unleSS a gross error was 
maqde . 

im anelivyeis Of the penetration problem only meen. 
Semeenergy was considered, with the sediment and object 
Meng a conservative system in the gravitational field. 
Meermal, chemical, electrical and atomic energies were 
mee rected. 

eee eciiod of Selucion Used with a Computer 

The numerical aleratiom Scheme wsed fo Solve 
Memetration problem was the same for each object shape. 
eee, LE divided the object into Sections or slices in 
Ene Same plane as the sediment Surface. Then the work 
Me@femodry LO Submerge each section and the work done by 
meen Of the submerged sections, the forces due to the 
Rieter OL added mass and displaced mass of the sediment, 
M@emoctactic Strength of the sediment, and drag in moving an 
@eeeemental distance equal to a section thickness were 
Setculated and subtracted from the kinetic energy and change 
mampotential energy of the object. Penetration depth was 
determined when the total kinetic energy of the object 
Semialeqd Zero. 

A description in mathematical terms follows and a 


block outline of the computer solution is shown in Figure 5. 
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Constants and mecessary approximations Fer seactuerieee naa. 
listed in Table I. The computer programs are included as 
Proencdixes C through H. 

Refer ste Pigure 3 for anxwexplanatron eimene 


See inicions that follow. 


ne eee ee nd 
k J | 
q Vv, [ 
v i 4 
bs fs 
“2 c rm 
wolrrs sonra X Position ] 
urs AYN AN inca oes SANDS) A WS 
DNA NS oes cnneseneaeaicen ele | CBS 
SESEESeiae DRae seg 
Stairs A iin 
f seem a te 
ts ‘J 
. f 
‘h Va A 
SIN Y pee \aS SS, ZN SAX Position 2 
PASEO Seems mises | aS 
SEES os. Jee Bisiinnnnnay 
CeIn tt osama aa ee eererney eet 
Teheametsttncaeecircatatatsncetetereners” g Telelefentsteterseteteterenente 
Mea ctcceretnscie, FMR atcaratetnetetereteteteses 


Figure 3. DLFINITIONS 
V. “saveuwage velocity of the object as liseraname- 


- Giseranecex trom DpOSition O LO post erenme 
(Assumes a linear velocity gradient) 
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v ~ average velocity of the objeGe as aie ence. 
distance x from position 1 tee pect iene 


M — Mass otf the object 


mM) added mass at section 1, plus the mass of 
Gee placed ft luvd aes cecect ion. 

mM, = added mass at section 2, plus the macs ven 
dra olaeoee tid iae Secu, one 

Am m, = mM. 


U —- work done 


U, ~ work done against buoyancy forces as object 
Leaticteamrr On PpOsit lon | togeeeserem 2 


Hee WOEMmeCOne against Sedamenk cltabic Stree Gama 
Lie ebweGus stranGSils fron, Osteo wie 
Peet Lon 2 


Us ~ WORKMCOne against drag forces vas the objece 
PhamsoltiSe: bom POStCTOn 1. bOmsOst monn 


OU the change in energy of the body as it UCeansi ce 
from POsSteicn 1] to position 2eenok anciweann 
ehange in potential energy 


Ey Be Otel emengy OL Ene CD CCl al geecuy ton 
ES pe tOucal, energy of the obj ect aleeesi «tom. 
moe nw eecheange an potential energy of theyobjecceme — a 


distance x 


From equation (22) 


2 , 
aa eta vari” (26) 

wT 

ws 2 2 
Wie gb em V5 - ip va Cae 
me follows that 
= 1 aye eee iy + + + 
eo = 1/2 mM, V5 - 2 m, Vy anol 2 a V5 U,, Ue Ug 


(2/3) 
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NOTE: 1. Am’ was accelerated from a State GConadte ton eee oe 
veloeity to the velocity at the ena on tame weeeeie. 
OVE Whrehs aie wac COMPU Led: 


2. Am’ Was assumed to be Zero aftew meena soe cern ac 
penetrated to its maximum dimension Since 
separation of the sediment from ehe Ob7eeruwes 
ap wee soem s Olt. 

Mies ncremental distance, x, is fixed, and the varicus Messesc 


eee Known aS is the work necessary to move the object from 


Meee dO! One to position two. Velocity at position two is 


Mimemown, Therefore from the work energy relationship 
E, = E, + PE - AU = 1/2 Mv? (29) 
2 1 2 : 
Memoining equation (28) with (23) and solving for V5 gives 
i —— a Se 
x + + } a = aaa 
ee 20088 Ey 1/2 my Vy ee eeu aa 
2 : (viewer ore m, ) 


meee 1teration scheme, when v, equals Zero, Ene Sun soma 


PieteenentS, Xx, iS the penetration depth. 

tiecawork Cone agaimstl buoyancy LOrCesS 1s ageecacwe. 
fOoLWward Calculation that is discussed in elementary fluid 
mrmrmrcs CLeCxts. in the computer iteration scheme where the 
eeremons Of an Object are taken to be slices with vertical 
eon. becomes; 


Z 


ees 1/2 yeh x (for the section being submerged) 


b 


plus boy & a (for all sections below the one 
being submerged) 


The work done agenimse Gdrag forces was taens ee. oe] 


Ua =e 2 pe A ye x (ea8) 
weere y ~ bulk wet density of the sediment 
C4 _ ~ arag coefficient for the approximate object 


co 


Shape that is below the sediment surface 


oo. 





A ~ the maximum cross Sectional area of the cbject 
below the sediment surface and normal to the 
Vy CllOG iy oC bos 


Vv ~— average velocity over the distemee me zca 
assuming a linear VariaquLOnw ae locas 


The work done against the sedimenc Stalbicwewecneen 


was calculated from the ultimate bearing Capacity mequamnaa 


Ue a C, N. Ac x (32) 
where Cy ™ a constant depending on footing Shape 
N. MecdecOMorant depending on the type "ob seaimaoms 
ancGiepoCmiiie roti ie 
c ~ shear strength of the sediment 


Nien tChespenceraling object. had -a flat plane omer a 
Meee cne Sediment Such as a cylinder impacting on end, the 
Pieees Sectional area, A, of the flat area and the shear 
eeemgtna, cc, of the sediment at the penetration level CEE 
Mmemierarca were used in the equation above. In the case of 
objects with rounded or sloping sides, a judgment had to be 
fereemwar d Cross Sectional area and representative Shear 
Meeemotni. The area selected was the maximum cross Sectional 
meeoemormal to the velocity vector at the surface of the 
sediment, or the maximum cross sectional area of the object 
if the area had penetrated beneath the sediment surface. 
This method of approximating a sloping or curved surface by 
a projected horizontal area appears valid for blunt objects 
[Meyerhof, 1967}. Since the shear strength characteristics 
of deep-sea sediments in the upper layers appear to change 
radically with depth in the sediment column, a shear strength 


for equation (32) was determined by selecting an average 








@epth according to the zones of plastic failure imelieaea 
equations (16) and (17). This was selected as one-half of 
the penetration depth until the maximum cross sectional aréa 
ime, been reached, then it became the depth of penetratren 
Minus one half the penetration depth when the maximum cross 
sectional area waS coincident with the surface of the sedi- 


ment as illustrated in figure (4). 







7 —— cross sectional area, A 


- 








Lt ea eee eee Toe 4h ‘ 
6: TT res ‘ \t v . eae warenora Ar on ae otenie ene soe sr ws a 7= xr Pete oo 
“\\ ANY woe a Rem \ Uy a 5 i lf ‘ ok. j po Yo 


= depth determining the value of ¢ 


cross sectional orea, A 







fe 
a «sediment surface 
"  SAVRTSNAG 


depth determining the value of e¢ - 


one half the depth when the maximum cross sectional 


area was coincident with the sediment surface 


Figure 4. DETERMINING THE DEPTH OF "c" FOR NON-FLAT OBJECTS 








The method described was adequate for curved and 
@onical objects and approximates a representative Ssereng en 
fo the theoretical zones of plastic failuwe in the sedimencw 
Bimalternative and more refined method of (selecting A anda 
would have been to use the differences in the cross sectional 
meee trom S€Ction to section and apply the value of € at vthe 
Semen Of Gach section. This increased the complexity G£ the 
Meeration procedure, hence was attempted only when there 
Meme discontinuities in the penetrating object's shape such 
femeohit Occur with a truncated cone. 

mie aduc@meance if Cach Tteration was chosen Oniy ete 
Pomserve computer time and approximate the expected accuracy 
Of the field tests. Larger and smaller iteration intervals 
Meme tested for effects on the final depth and as long as 
the interval was less than one-half of the distance between 
shear Strength measurements, convergence was always apparent. 
@eeedea Linear shear strength profile, the interval is limited 
eimeveby the desired accuracy of the solution. 

Wieweclarrons £Or accdcd maSs are Given in Taotewa 
feo chart of the basic program iS given in Figure 5. 

Appendix A presents a set of curves which qualitatively 
Mmftenrcate the relative influence of each force in the pene- 
tration equation. They were obtained by using the sediment 
Pm@ear strength profile typical of the world's deep-ocean 


basins as determined by Hoag [1970] and shown in Figure l3. 
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Made of motion 
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Radlus Moving 4 2 
m, = k a 
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Planetary ellipsoid 


or oblote spheroid 

: = eine ee 
moving poroliel to Ui 4 e~(sin lay ae 
Its polar oxis yi a i wee a 





Etliptic disk meving 
perpendicular 10 its 
plane 


m, =ha ( 3 pmo b) where 
ake 


k4 = (07 oe 
La ine a2 . 
5 [e sin?@ + cas*6) dG 
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Table Il. FORMULAS FOR CALCULATING ADDED MASS. 4 


tHoerner , >. Fo, Piula Dynamic Drag, Practical intouiaemoume 
Pemeocynamic Drag and Hydrodynamic Resistance, Chap. If, Pub- 
Pished by the author, Midland Park, New Jersey, 1958. 


STerzachi., Karl and Peck, Ralph B., Soil Mechanics [iE iigme ae 
meee actice, p. 1/72, John Wiley and Sons, 1948. 


Seeverhot, Ge Gea, ‘The Ulcimate Bearing Capacity of VWeage 
meeeeo Foundations’, Proceedings, Sth International Conter— 


eee Oil Mechanics and Foundation Engineering, v. 2, p. i0o, 
Paris, lL96l. 


* saunders, Harold &., Hydrodynamics in Ship Design te 
meee?! , The Society of Naval Architects and Marine Engineers: 
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Figure 5, OUTLINE OF COMPUTER PROGRAM 
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GC. MODIFICATION OF THE METHOD FOR USHA ier. 

The method was ultimately developed for general use at 
feemewith aS few restrictions as possible placed on si 2enon 
Shape of the object and on the type of sediment. However 
Peeeonly data avaibable to test and verify the methodmwas 
meet alr-to=-simulated deep-sea sediment case. 

Here are thes relatively simple alterations thatweeneule 
memmacde to the computer program before it can be used in a 
meeeer-cLo-sediment case. The three alterations should act to 
G@em@ease penetration depths at higher velocities and their 
Meee! effect is to increase the importance of an accurate 
Pee-me@ination of the shear strength profile. 

bemect, the added mass effects will change by the ratio jor 
Mmemeensity Of sea water to the density of the Sediment. This 
Bamecactiy be accomplished by dropping the unit value of Sea 
memeer from the bulk wet density term in the programs. 
Additionally, the added mass of the object free-falling 
through the water inust be added to the total kinetic energy 
Mimene OOJect at impact and a method for calculating the 
Seeeeets OL this entrained water mass taken into account. 

Beecond, the buoyancy term will also change by the ratio 
femme GensSity of the sea Water to the density of the 
Becaiment. 

Mmrrd, the drag term will change by the ratio of Sea 
water density to sediment density. Also a drag term must be 
Mm@eluced for the part of the body that projects above the 


ieee : 
sea-sediment interface. The easiest way to accomplish this 








would be to include a drag term for the whole object and add 
peorag term for the portion of the objec@isubmergecms een 
Peqaiment. 

The three modifications will serve to reduce the influence 
of added mass, buoyancy and drag and cause the sediment 
mearerc Strength term to become even more Gominant than in 
the air-to-sediment case. Keller [1964], used the standard 
Meemerng Capacity equation, only, to predict the initial 
Mmemeereoment depths of a rectangle, Square, and cylinder into 
soft bay sediments. The results were poor and he concluded 
meee OoLObLems with determining an accurate strength of the 
Pearment were the reason his calculations always underesti- 
merece the settlement depth. 

To correctly use the method at sea, several cores should 
be taken and their physical and mechanical properties averaged. 
Peeraging Sshovld reduce errors caused by sampling and testing 
Miperutoances. The effect of Sensitivity, aS it applies to 
the shear strength used in computing penetration depth, 1s 
Mmoecertain. Possibly the remolded shear strength iS a better 
Malue if it is assumed that impact shock effectively remolds 
~he sediment to an unknown extent. In sensitive sediment, a 
range of predicted penetration depths should be calculated 
meeorgding to the degree of sensitivity. 

Mme orientation that the object haS aS it fixSt Starrs ee 
Menetrate has a significant effect on total penetration. 
After free-fall in several thousand feet of water, the 


object may still be tumbling or oscillating if it is unstable. 
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menmida [1966] found that a Significant portion of a falling 
body's energy may be in its angular momentum which is not 
included in this method. 

For these reasons when predicting penetration depths, a 
memputer routine Should be used to solve for a range of impact 
Merocities, shear strength profiles, and object orientations. 
The calculations can be done by manually using larger inter- 
ValS, aS was done to check the first computer solutions, but 
Mee energies are being equated and the method is iterative, 
extreme caution must be used Since an anSwer, correct or not, 


meways results. 
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IVs RESULTS WO EAPERI MER is 


A. THE PROBLEM OF MODELING IN SOILS 

Some problems with testing and modeling the dynamic 
behavior of soils are discussed before reviewing the data 
Mmmeaole to test the penetration prediction method and the 
Mesults of the tests. 

mee field of soil mechanics has traditionally been con- 
Meemecd With large scale construction operations. Many of 
meme Lnods used to calculate the stability of slopes and 
Meroe tion designs are based on ultimate collapse criteria 
and much of what has been learned about soil behavior was 
Memavea from studics of catastrophic large scale failures. 
meemconstruction projects are seldom undertaken without 
Memeensive, full-scale, field tests of critical portions of 
Peeeeeoundations, and mocern buildings are carefully instru- 
eee tO Monitor their settlement and reaction to seismic 
Beeeturbances. A good indication of the state-of-the-art of 
Seeeeechanics is implied in the large safety factors 
meeemmended in texts and building codes. The lack of under- 
memaing 1S, in large measure, due to the inability to 
Parameterize and model a substance as complex as a soil. 

When the problem is narrowed to the field of saturated 
deep-sea sediments, the problems of modeling are somewhat 
meaqucedq. A modeling approach to investigate the dynamic 
response of soils to impulsive loads was conducted by Westine 


me766| who found that a satisfactory case could be made for 


eel 








a model with saturated clays, but granular materials with 
epgles of internal friction were impossible to Simulate, 
Since dilation and pore water pressures were hard to model. 
Mmeolich [1967] attempted to model the dynamic behavior of 
Submarine sediments, but unfortunately he passed away before 
Sout shing his final results. 

mmere are several areas where modeling dynamic pene-— 
Mmeation phenomena in deep-sea sediments may prove to be 
mmoOossible. 

mercer dynamic loads a saturated deep-sea sediment may 
Penave as a bingham body. The inertial resistance of the 
Sediment to deformation coupled with the threshhold stress, 
Meotatic resistance of the sediment, may be very hard to 
Memarate if the mass anc therefore, inertia, of the pene- 
eating body is small or on the same order of magnitude as 
Memo tatic resistance and inertial resistance of the soil. 
Meena logous problem, but Simpler by one factor, is unsteady 
muta flow at low Reynolds numbers where inertial and viscous 
Mmeeces are the same order of magnitude. In the dynamic 
meee Of Soils this problem can be circumvented by full- 
Scale tests using objects with large masses. 

A Similar problem involves surface effects and the 
MmeectOn of deep-sea sediments to an object. The effects 
ee adhesion have been shown to be influenced by polymer 
Meetings, Surface roughness, surface chemical composition, 
Becetera, by Erchul [1969] and other investigators studying 


Mme driving on land. Surface effects may control penetration 
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Semi ficantiy af the object has a largevsuLrlace sarecamim 
Belation to its total mass. A means Of minimizing adhesion 
Meeects would be full scale tests with large masses. 

Ef a study's purpose is to learn more about adhesion, or 
mass effects, then modeling may be the best approach. However 
MeecapOlation of modeled test results to objects larger by an 
Beaer Of magnitude may be invalid. 

manaliy, when attempting to model a sediment composed of 
Meey particles which are very small, to reduce them further 
Meee attempt at modeling the interparticle attractive forces 
and pore water viscous effects, may be impossible. Also, the 
instruments available to measure the modeled sediment strength 
fmmeouw Nave the required accuracy or precision for a reduced 


scale. 


fee DATA AVAILABLE 

The Naval Ordnance Test Station, Pasadena, California, 
eee Led a large pit from sheet piling, 37.5 feet in 
Giameter, and filled it to a Gepth of 16 feet with a very 
Mek sediment from Seal Beach Lagoon, California. The pit 
@eertcead to conduct tests on the holding strength of anchors, 
Meeakout forces of objects imbedded in the mud, and experi- 
Mecs With dynamic penetration. The penetration experiments 
Mumma rica out under the direction of Dr. R. J. Smith in 
September, 1965. The penetration experiments are described 
below. 

Meee sctiption of Obyects and Tests 

Figures 6 through 11 show the six different objects 


that were dropped into the test pit. All six were fabricated 
oe 





The Spheve wes Madewo ane 
heavy "seeeT plete ram.s 
WElLGNeG S560 pound cman cede 

A padeye preovidedema cero en — 
MeCnE POI e: licmeninen 





Figure 6. SPHERE 








The horizontal cylinder was made from a 55 gallon 
Pr drum, itt weighed 540 pounds in azz, it was 
meeced bya two penmmant bridle wren one-pennane 
Pexded at each Lip. 


Figure 7, HORIZONTAL CYLINDER 
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Figure 8, 


The cube was fabricated from 
feet Plate and weighed 702 
Meunds in air. It was lifted 
mea LOUr pennant bridle to 
mel CcOrner of one face. 


Figure 9, 
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The cylinder-on-end was also 
made froma omc Vile neane: 
drum. For the first tests 
it weighed 540 pounds in air. 
hater a@iiac, sounder te 
was WEClCeGG 367 or VEie simaom@dn 
Che GOruM ae eene onan ema 
DeEReEEL ale maki NG Vic soe wid 
LO Leduce —Ehe amount: Gfeagur 
Biat Might Nay eowoeen teeppee 
there. The plate increased 
the Weight to So) peundee rn 
air. ~The drum was Iie diaiy, 
a three pennant bridle attached 
iOmnineo | emia r cid, 


CYLINDER-ON-END 








CUBE 








= = is rth, 
Se The rectangular para- 


Piles peace) Cincom 
ny ip Ea @unaenon li ane 
Ti. Was Ons ruc recom 
steel plate and i197tedq 
ite with a four pennant 
. oe DELLE” CO" Our Comme ms 
Al 3 Oc omen e - Sider 


Figure 10. RECTANGULAR PARALLELEPIPED 


mee COne was fabricated of 
Meeer Plate and weighed 
mee, POUNdS in air. It 
Meer ai ted by a three 

meet pennant to the flat 
meade and dropped point 
down. 





Figure 11. CONE 
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mith a steel shell for Structural FigiGi a sad ee gee. eames 
meme shell was completcly filled with £uesh water for weigqic. 
Mne weight could be chanced if desired, however for pene— 
meation tests they were pressed full at all times. Dimensions 
Bf the objects are shown in inches. 

The objects@wexrc suspended sby a Crane ovemweReselectce 
Seeerons in the mud pit. They were released by a quick 
memease hook that was tripped manually by pulling a line from 
mnie Side of the pit. The objectS were released at different 
Meepents above the mud-air interface to allow them to reach 
mmmeetics of 5, 10, 15, and 20 £fcet per second as the 
object's lowest portion contacted the sediment-air interface. 
mmeezcrO velocity test was conSidered to be a zero release 
Mme, and the objcct was released with its lowcst portion 
Just in contact with the sediment-air interface. 

Slow motion pictures, 32 frames per second , Gf vainos. 
feeeeee penetrations were available. They were carefully 
memcda1cd in a frame-by-frame projector to double check the 
Seeecata Sheets for height of drop and penctration depth. 
meeemovies also gave an indication of the Stability of the 
Memeece aS it penctrated. Remarks about each penetration are 
Men in Table X. Times for penetration and times for 
meeound weorc estimated by counting the frames of the movic 
film for cach event. The movies confirmed the observations 
meechiuang [1970] that air trapped below a body impacting on 
Peeuid does not alter the penetration depth. Any air that 


Meet rapped bclow the body uSually vented by bubbling up 


oy 








meanaside the object after it had come (uOsemeee eon): win 
one or two caseS was it observed to form a jet with a 
Beenificant amount of entrained mud. 

The penetration depth was determined by ie 
observation of marks painted on the objectS which were scales 
Seeded into one quarter inch increments. The depth refer- 
ence was the undisturbed surface of the pit. Depth obser- 
Memerous were made at each corner of the rectangular objects 
Meeeat a diameter of the axisymmetrical objects. The pene- 
meation depth was considered to be the numerical mean of the 
Beeral measurements. Any exceptions to this method and 
Meeres Of tilt at the object's final rest position are noted 
emma remarks column of Table X. The exceptions occurred 
een the cylinders-on-end tipped completely over after coming 
Memrest. The penetration before tipping was used if it 
meme pe read from the movies with an accuracy of one half 
mnch . 

meee ccription of the Sediment 

mere was a real problem finding a Suitable material 
Memsimulate a true deep-sea sediment, compounded by the 700 
Mere yards of material required to fill the test pit. The 
Meeerial selected was Lagoonal mud, described as “unconsoli- 
Mmeed grey sediment with medium sand sized mica flakes", from 
the tidal flats of Seal Beach Lagoon. Tables III through Vv 
Contain detailed analysis of its mechanical properties, and 
mineralogical and chemical analyses of the lagoonal mud from 
Mees taken from the borrow pit. A year after the pene- 


tration tests were completed, Muga [1966], while measuring 
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TABLE Il. MECHANICAL PROPERTIES OF THE SEDIMENT IN THE TEST PIT™ 


Depth (in) 0-3 4-7 
Color (GSA No.) 5GY 6/1 5GY 6/1 
Odor H,> H9 
Bulk Wet Density (gm/cc) [e452 1.559 
Vane Shear Strength (psi) 0.364 1.495 
Remolded Strength (psi) 0.364 1.736 
Sensitivity 1.0 O27 
Water Content 82.9 66.3 
Specific Gravity of Solids 2.914 3.07 
Dry Density (gm/cc) 0.794 Cayo? 
Void Ratio 22676 a2 
Porosity (%) 72.8 69 .0 
Saturated Void Ratio 2.416 2,002 
Organic Carbon (%) loco roe 
Percent Sand 6.0 Za 
Percent Silt 30.0 S,0 
Percent Clay 64.0 80.0 
Median Diameter (mm) 0.0020 0.0018 


Note: Samples were obtained with a two and one-half inch core liner from 
the site where the sediment originated. 


Beamples ColleeGeed= by k. J. Sten and analyzed at Une 
Seer CUVil Emqineering Laboratory in Port Hueneme. 
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TABLE IV. MINERALOGICAL ANALYSIS OF THE TEST PIT SEDIMENT® 


COARSE SILT TO SAND PERCENTAGE 
Quartz WS, 39 
K-Feldspar PeeC? 
Calcite OG 
Opal 1.47 
Muscovite (illite) 11.54 


Total = 39.45 percent 


ACTIVE CLAYS 
Kaolinite 16.3] 
Mixed Layer Be oe 








Total = 35.58 percent 


INACTIVE SILT TO CLAY SIZES 


Gosthite CF 316 
Halite 4.53 
Organic 4,4] 


Total = 14.77 percent 


Note: The mixed loyer was composed of 67 percent montmorillonite and 33 
percent chlorite. The muscovite was chiefly terrigenous (clastic) and occurred 
in coarse silt~sized grains. The analyses were made by a combination of x-ray 
diffraction, chemical, and petrographic methods. Infrared absorption analyses 
and differential thermal analyses were used to obtain supplemental data. The 
values shown are recomputed to 100 percent and are reliable te about +15 
percent. 


Sinelan, J. Aj, An Report "on VUineralocy Cimee-an Peecon 
Beaiments, p. 10 unpublished report to the Naval Civil 
Mmegineering Laboratory, September, 1967. 
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TABLE V. CHEMICAL ANALYSIS OF TEST PIT SEDIMENT” 








OXIDE PERCENTAGE 
Si 0, A8 34 
Al. 0, iomZo 
Fe,0,, 7 88 
MgO 27 
CaO 2 
KO 2.74 
Na,O 3.02 
HO Zhe 
Lae Tay, 
Total a 
Surface HO 73.01 


‘Whelan, io © Neer oa 


Hol 





TABLE VI, RANGE OF VARIATION OF TEST PIT SEDIMENT™ 


"The main feature of the soil is its highly variable properties, being 
both nonhomogeneous and nonisotrophic. Some of the variations are indicated 
below. 


RANGE OF VARIATION (42 TESTS) 


MINIMUM MAXIMUM AVERAGE 


Vane Shear Strength (psi) 0.0 ree a0) 0.530 
Remolded Strength (psi) 0.0 0,861 Oneod 
Original Water Content (%) AS2o oe 6345 
Void Ratio 1.104 2.126 le71G 
Porosity 5255 68 .0 Ole 
Liquid Limit 41.0 eo: 0 56. 1 
Plastic Limit 77 6o 43.4 34.2 
Liquidity Index 70 Zoo 136 
Percent Sana 1.7 Lae 1052 
Percent Silt 30.8 94.4 40.5 
Percent Clay 35.8 66.0 46.4 


"In addition, soil samples were subjected to numerous triaxtal, permeability 
ond viscosity tests. From these tests it was determined that; (i) the yield 
strength of the ‘soil’ is approximately 0.5 psi; (ii) permeability is practically 
nil; (iii) cohesion, c, is almost zero; (iv) angle of internal friction is almost 
zero and (v) viscosity is independent of the rates of shear or deformation 


fe. tpm)." 


ease , feo pe erCan Out erOrees, "O.64,5U..0o. Nana lysGa al 
Emogineering Laboratory Technical Note N-863, Port 
Beememe, September, 1Y66. 
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break out forces in the test pit, conducted more tests to 
Meaty the mechanical properties of the Sediment] Table 
mis a Summary of Muga’s tests. It iS Obvious that the 
memeonal mud was not uniform but characterized by variability, 
but a major cause of the variability was undoubtedly the 
Meet OF Sampling disturbances and continual reworking. 
Meee alSo stated that the mean viscosity was 300 square feet 
Der second. This value represented a series of meaSurements 
Meer With a commercial rotating-type viscometer. He did not 
mttcempt to correlate the viscosity with the vane shear 
Meeeength because of the "marked variability in the soil 
Meee rties’. Since the sediment contained approximately 

40 percent silt to Sand sizes, there was certainly some 
emer OF internal friction caused by these large sizes and 
Semeeewas probably the reason Muga found a high viscosity 
Meeeourement. It alSo appeared that there were more coarse 
Mein Sizes present a year after penetration testsS were con- 
SBucted which indicated possible contamination with other 
eee tial. For this reason Muga's values were not taken as 
Mummece |yYy correct, but they were included to illustrate 
Sem—~evariability of the sediment used for the tests. 

Figure 12 shows the mean shear strength profile of 
mae sediment from 550 meaSurements. Although the mean 
Merength profile is not representative of the average shear 
Berencth profile for the world's deep-ocean basins shown 
Mmeea SOlid line, Figure 12, from Hoag [1970], it 1s 


remarkably close to the shear strength profiles of the edges 
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DEPTH (inches) 


SHEAR STRENGTH (psi) 


| 
18h 
24 
30 
36 
42 
48 





Are mean 
ee one standard deviation from mean 


-average shear strength from 
14 cores from water depths 
greater than 6,000 
[Hoag, 1970] 





Figure 12, TEST PIT MEAN SHEAR STRENGTHS 





Of the Zigsbee Abyssal Plain as shown by Bryant and others 
mro67). 

The numbering SyStem designating Lhewleeae renee. 
Peation, where each penetration test was conducted is shown 
Semrigure 13. The shear strengths at each Station were taken 
Mmeang a two inch by five inch vane. The vane was operated by 
See cLrom a punt floating on the surface. It was driven into 
feemeoediment in Six inch increments and torque applied by 
meano With a standard torque wrench and the maximum torque, 
moenverted to shear strength, recorded. Not all of the vane 
@eeees Were conducted in the presence of an engineer and some 
Meee done by different personnel. Generally the shear 
Memeeigch profiles are consistently Similar, with some erratic 
Mees, but there is some correlation between very erratic 
Meeles and the date of tests, which indicated operator 
feeeeeand care may have had an influence. The method of 
Seeing Shear strength values greatly influences the value 
Mmeained according to Terzaghi [1960] and it was difficult to 
Mee consistent readings from an unstable punt with no control 
Meee Vane's rotational speed. 

Other methods of determining shear strength profiles 
Mmoetne mud pit were considered, but rejected for the follow- 
Mme reasons. In 1965, the existing penetrometer probes 
related resistance encountered to shear Strencth by ciipleicam 
Beables or graphs, and the very low strength of the mud pit 
Was out of the empirical range. The use of these devices 


would have reguired modifications of the instruments and 














Siw? “aa Cee, <n” sui eee? ana 


DEPTH OF PIT = 16! 


NOIE> Small cirete indicates tie boecation of 
Sub-lettered clear strength data pointe 
around a maine station. 


Figure 13. LOCATION OF DROP STATIONS IN THE TEST Pit 
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eerensive testing to obtain empirical relaciveonsi1 3s. ooecenee 
mae sediment in the mud pit probably had some Slight but 
mesnown angle of internal friction and the vane test mea-— 
Pures a cohesive shear strength which is influenced by the 
Pere Of internal £riction in a minor ways However, the 

meme test seemed the most practical test for shear SCECnG ens 
_umeeres and in-Situ testing. Therefore, vane shear strengths 
Seeeesced exclusively to keep the resultS comparable to a real 
Situation. 

The horizontal variability of the observed shear 
Mereengths at six inch depths is shown ina plan view of the 
mead pit, Fiyure 14. The measurements taken on different 
Gates in the same location also show a spread in values. 

ion Vabidellity OF Liesshear SuLengii ell. uneumia mene 
Serer eLObably not random for the following reasons. 

i The huge volume of material required tc £111 che 
Meeenecessitated handling by large construction equipment 
Meee COUlLd have scooped up isolated chunks of material and 
Memmeved them to be placed in the pit essentially undisturbed, 
Matle most of the sediment was greatly disturbed during 
Meeevation, hauling and placement. The very high shear 
Merengths have been attributed to the relatively small vane 
moming in contact with one of these hard chunks. The chunk 
merc Mot necessarily affect penetration depth unless an 
Bojyect was dropped very close to one because the sediment 
Meomancompressible and moved as a unit. Most of the objects 
emeeelarge compared to the chunks or vertical distance 


between anomalously high observed shear strengths. 
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NOTE. Mores than One value Indicates measurements 
were taken on different dates. 


Figure 14, SHEAR STRENGTHS AT SIX-INCH DEPTHS PLOTTED 
ON A PLAN VIEW OF THE TEST PI (psi) 
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2. Althougn if was not noted I Baie ona) 21. io. 
mae sediment, some organic material or fe@geign objects may 
mave been present Since no Special precaugions were taken 
Memeoteve the sediment while it was being placed. 

3. The extremely low walues of observed shear 
Meeength may be attributed to the method of handling the 
Beemment after a penetration trial. The holes made by the 
objects and other experiments with the sediment were smoothed 
over, filling depressions with a shovel and rake. Then shear 
Mmeemoths were taken at each station and penetration drops 
memoucted at a later time. After all the stations had been 
meed, the procedure was repeated. Since other tests were 
meee COnducted concurrently, the hole filling and smoothing 
Seed CONntinuous operation. Undoubtedly some very soft 
material and water accumulated ab ENC DOtEonS Lr chechole = 
Mee cits accounts for the anomalously low shear strengths at 
Meioiderable distances below the surface. 

mie low sensitivity of the sediment before Sampling 
Seeecne COntinual reworking it underwent while tests were in 
Meee ss, led to the conclusion that the sediment used for 
Muece sts had a sensitivity of one and that the observed 
Meeettlu shear strength of the sediment was a good repre- 
Meteo ciOon of its actual in-Situ strength, if the problems of 
Meee Operator error and the anomalous points due to hard 
meunks or soft spots could be resolved. 

Memetiinately, enough data, 716-poines =exloeeos FO Man. 


Mmeeeatistical analysis and a méan strength profile of the 
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mea pit. Table VII includes a matrix CE Valier yeseeriee 
Shear strengths. The anomalously high and low values were 
Meyentified by eye. Since the profile should have been 
mearily smooth with all the reworking the sediment underwent 
Sing placement and tests, the hard and soft Spots were 
feeeey Selected by visual inspection. The values considered 
anomalous are underlined in Table VII. 

mre mean and Standard deviation Of vsnean Strona tomas 
Meeeinction of depth, both with and without the anomalous 
Values, are listed in Table VIII. It should be noted that 
@eenwanomalous values often cause the standard deviation to 
Meee than double. The mean shear strength profile of the 
moe pit and one standard deviation calculated without the 
mromalous points are shown in Figure 12. 

mercaict Station tie Beonalore points were adjustoed to 
Bec a Smooth profile that generally was within one standard 
Seeeration of the mean. The adjustment process was again 
Meojyective, as the precision of the original data did not 
Smeranit a more sophisticated approach. The observed shear 
Meeenoths for each station were plotted with the mean pro- 
Pele. Then the anomalous points were moved to fit a general 
BmoOth trend of points within one standard deviation from the 
meee «6 The figures in Appendix B illustrate the adjustment at 
Macn Station and the solid line is the actual shear strength 
meorile used by the computer solution to predict penetration 
mepth. 

The shear strength and standard deviation value at 


the air-sediment interface was taken to be zero, in that sea 
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TABLE VIII. MEAN SHEAR STRENGTHS AND STANDARD DEVIATIONS OF 


ES Tate 
DEPTH MEAN SHEAR STRENGTHS ONE STANDARD DEVIATION 
(in) (psi (psi) 

6" es | er) ee .38 
2" 73 (60) 79/73) eye 4 
13" SEL ates) 47 . 84 
24" U0) ey 63 le? | 
30" ee) 2.22 5 62 eve 
36" pes oo Z a -/Q 1.94 
42" Vea) S78) 60 ter 
48" BENE eal 74 1.60 
54" 73 (95) ee evs oi 7 
60" 3) oe J | 1.8] 
66" 2.92 98) anh) oe Za 
u2" 2092 04) SO ev Or 1.49 


Note: The first columns under the shear strengths and standard deviations are 
values computed with the anomalous points discarded. The second columns are 
values computed with all the data points. The subscripts in parenthesis are the 
number of data points. 
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Mecer waS continually added to the test pit to offset 
Meemooration losses and the top of the pil wasseentinuall, 
meng reworked by raking. In general, a@ Virgin UnarsTturoed 


meeement Should have a finite strength at the suriface. 


fen RESULTS 

meble IX summarizes the observed penetration depths, the 
Meeemceted depths, and tabulates the difference between ob- 
Mmewed and predicted. A numerical analysis of the differences 
Meeeween Observed and predicted depths was not made because, 
mumeecew exceptions, the method developed accurately pre- 
Gicted the penetration depth within the experimental uncer- 
mm—mmetes Of depth measurement and shear strength accuracy. 
meermethod overestimated penetration 17 times and under- 
Seer ed penetration 22 times. The largest differences were 
Mememestimates of penetration. In all cases of differences 
Mueeteer than four inches, the prediction had to rely on an 
Memustcd shear strength. This indicated a problem with 
Meeesuring the shear strength at that station and implied the 
Meeeeence Of a soft pocket in the sediment. The large dis- 
MeeDancies also occurred with impact velocities of five feet 
mer second. 

Meme X sets forth the penetration data along with notes 
eee cOm movies of the tests. The time-~to-~penetrate was 
Bbserved by counting the frames of the movie film from object 
ampact until all downward motion had ceased. The rebound 
time was also observed by counting movie frames. The rebound 


distance was estimated from the pictures and is believed 
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Table IX. SUMMARY OF RESULTS 
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Oy ea | PENETRATION 
VELOCITY OBSERVED PREDICTED DIFFERENCE 
(ft/sec) (in) (in) (in) 
BPRERE (wt. = 561 Ibs) 
0.0 eh 1Oe7s fe 
0 eee le =<4RS" 
10.0 14.0" 14,0" o207 
Hess. 24.2" 2A FEA) oh 
20.0 AO 20r oe TOs 
ee ome iess ~4,5" 
Ho). O esi" 278 ie =O02 
ie, O 920 Zoa0 14 Ue 
zy. 0 le vom 20.4" eas 
HORIZONTAL CYLINDER (wt. = 540 Ibs) 
2.0 D200 ole + 0.5" 
> 0 tee: Ton - 4,0" 
HO. O rot loo oa 
io, U eos AAO e O20" 
20.0 lose leno: - 0.]" 
20.0 las 14,7" Cz0! 
CYLINDER ON FLAT END (wt. = 540 Ibs) 
>. 0 14,0" Oe - 5,3" 
ie), 0 oreo PA mae a2.7.' 
ES, O Ise" Ia - 0.8" 
20.0 tee V6.2 me ee 
0.0 8.0" Sa4 +0,4" 
S, 0 4.2" A,4" ae = 
10.0 eS ee Ne” - 2,0" 
ie. U 14.0" 257: - 1.3" 
MoO 18.0" IA Or ~ 3A" 
20.0 OAs KOvou = 262" 


TABLE IX. (CONTINUED) SUMMARY OF RESULTS 


WBIECT PENETRATION | 
MELOCITY OBSERVED PREDICT iE DIFFERENCE 
(ft/sec) (in) (in) (in) 
CYLINDER ON FLAT END (wt. = 561 Ibs) 
5.0 0 een 12.68 ee 
10.0 14.0" 14.0" 0.0" 
ie:, O Sees” 20.4" lca 
0.0 nee onlin Ole 


SUBE (wt. = 702 lbs) 


0.0 4,5" Bf Ono! 
a 0 i 4,4" -1.3" 
10.0 So" es ota 
15.0 Von i268" =n 
2050 oe 14.9" =/2 61" 
RECTANGLE (parallelopiped, wt. = 1014 Ibs) 

0.0 S50" Sue OR. 
>, o.0 7a Oe 
10.0 wale 94" Ore 
15.0 14.8" 14.7" anne 
15.0 loon 10! One 
20.0 PANTS Zale We 
CONE (wt. = 1127 Ibs) 

UG S50 | 50.0" -5,0" 

5.0 Jo SiOy =2.)" 
10.0 49 5" 54.0" +4,5" 
15.0 | 65.04 64.0" Or 
20.0 O2.On Do aee =o. on 


ow 





Peeurate to one inch, Penetration depth were: {eee ene 
Siepect's final resting position and did not exclude the 
mebound distance. The percent of ultimate bearing capacity 
meme alculated as; the ultimate bearing Capaciey of tEhe 
Beereaent at the depth the object came to rest, divided by 
Memeactual force applied to the sediment by the object at 
Semeermal resting poSition. The shear Strength adjusted 
meu indicates whether or not an anomalous shear strength 
menue was adjusted which directly affected the calculations 
Smee netration prediction. 

Mhe time-to-penetrate measurement waS done in an inexact 
Manner, but is remarkably constant when the range of 
Mebocities and different object shapes are considered. Its 
Meet ive constancy in relation to the other parameters may 
indicate oct tc lomrelated tO ad Unley oases peOper ey som 
Meeesecaiment. Schmid [1966] observed the same phenomena 
mee cesting small penetrometers with different impact 
Meerocities. 

Mire cebound time and distance were very difficult to 
Meetmiate due to the short duration and fact that splatter 
Bemetimes partially obscured the object. Some of the rebound 
Of the object, and apparent elastic behavior of the sediment 
mee Originally thought to be caused by a buoyancy effect. 
However the cavity formed by penetration did not collapse 
Meer! rebound was over and the object was essentially at 
Mest. Subjectively it appears that rebound is a function of 


Memelastic property of the Sediment. 
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The percent of ultimate bearing capacity that was reached 
Mecne different objects seems to relate fue impeces clocme 
mmol the penetration tests the objects were lett in plage 
hor at least 24 hours and in some cases up toa@ght days and 
mm@en checked to see if they had settled further into the 
meement. No further settlement was observed in any case. 
Meee impact velocity of five feet-per-second the objects 
Memoistantily penetrated to a level where 300 percent of the 
Meeeimate bearing capacity was available to support the object. 
Seeoeimolies that a structure being emplaced by a free-fall 
Meemnigque would only require a velocity of five feet-per- 
mone CO generate the safety factor of three recommended in 
Memeeemechanics in a Situation where the ultimate bearing 
Meoacity equation is used. 

memmeoserved in the movies, any air entrapped under the 
Meee nad negligible effect on either rebound or pene- 
Meeeton depth. Rounded objects such as the sphere and 
MemezOntal cylinder did not behave differently from the flat 
Memeo OOj;ects tnat did entrap Some air. 

Mieeeering and pile up of the sediment was slight. fin all 
Meoes, even with the rectangular parallelepiped impacting at 
memreet per second, the crater lip was less than 2 inches 
migh. The sediment moved laterally away from the objects and 


mOC up. 
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IV. CONCLUSTONS AND RECOMMENDATIONS 


Pm rational, phenomenological, theory Haseeeenmapolieadmas 
mee problems of free-fall penetration into a Simulated deep- 
Sea sediment. The theory correlated well with full scale 
Mots. A more complete solution will probably only be 
ert aole after more is understood about the mechanical 


Beeoperties and dynamic behavior ‘of deep-sea sediment. 


a) 6 CONCLUSIONS 

Bemememethod predicts the penetration Of f£ree—tai ling 
Objects into weak, saturated sediments for an air-to-sediment 
case within the accuracy of the present state-of-the-art 
meeeeraquies for measuring the sediment mechanical properties. 
Meesiiould work equally weil in a water-to-sediment case with 
meme modifications. 

fee impact duration time, from the instant of contact 
Merl the object comes to rest, was observed to be rela- 
Bavely constant and may bc a unique property of the dynamic 
Mem@evior of a sediment type. The duration time was rela- 
Bavely constant even though the penetration distance changed, 
Mupemmeoj;ect Shapes were different, and the kinetic energy of 
Meemecatling objects ranged from O to 5,000 foot-pounds. 

The method may be used with larger, heavier objects but 
meution must be used in extending the method for use with 
Meet ler objects since viscous drag and surface adhesion 


Merects may be important for smaller masses. 


oy 











At high impact velocities the effects of addeamiac- ean 
mem drag cannot be neglected in computing penetration depth. 
Mae accuracy of the constants and variables in the standard 
Meeeing Capacity equation have the most Control over pre— 


meted penetration depths. 


B. RECOMMENDATIONS 

me method shemld be experimentally verified by full 
Mme tests at sea or by controlled tests in a large tank 
Mmemore USing a uniform sediment whose mechanical properties 
Meenrknown. Experiments of the scale needed would require a 
meepor Cffort but could include investigations into other 
mmenown areaS such as the free-fall attitude and stability 
Memevarious shapes in water, the verification of terminal 
Merocitics of free-falling unstabilized objects in a water 
Meeeumn and the study of breakout forces and salvage methods 
Meed to retrieve the objects. 

Pemenevyer an object is recovered from the deep-ocean floor 
Mey Cffort should be made to obtain photographs of the 
meer. The sediment around the object Should be extensively 
memo lLead and completely analyzed to allow verification of this 


Meeod Of predicting penetration depth, 








APPENDIX A 


MORCES ON OBJECTS PENETRATING A TYPICAi weep ivan. 

Byopendix A includes a set of Curves tihici qualitative, 
Mmiictrate the relative influence of each force in the 
Benetration equation. They were obtained by using the 
meaqiment shear strength profile typical of the world's 
deep-ocean basins as determined by Hoag [1970] and shown 
Semeegure 13. The curves are normalized for depth of pene- 
Bration on the abcissa and the maximum force encountered on 
Meeocdginate. Time is independent of the calculations and 
Mupeeetape of the curve near the origin is largely a function 


See the iteration interval. 
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me- Shear Strength prezives apweachiee veo (stati ea ore 
Berustrated together with the adjusted profile used in the 
Meeeration prediction computation. The profiles are in the 


ame order as listed in Table VII. 
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